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We formulate an effective model for B-B' site ordering in double perovskite materials A2BB'06. 
Even within the simple framework of lattice-gas type models, we are able to address several exper- 
imentally observed issues including nonmonotonic dependence of the degree of order on annealing 
temperature, and the rapid decrease of order upon overdoping with either B or B' species. We 
also study ordering in the 'ternary' compounds A2BB'i_yB"^06. Although our emphasis is on the 
double perovskites, our results are easily generalizable to a wide variety of binary and ternary alloys. 



I. INTRODUCTION 

Double perovskite (DP) materials of the form 
A2BB'06 are being actively studied^ on account of their 
unusual electronic and magnetic properties. In particu- 
lar, some double perovskites, e.g^ Sr2FeMo06, show half- 
metallic behaviour, high ferromagnetic Tc, and large low 
field magnetoresistance. Remarkably, there are also dou- 
ble perovskites that are antiferromagnetic and insulating. 
The metallic and magnetic character depends primarily 
on the choice of B and B', but is also affected by the local 
ordering of these ions. 

While it is desirable to understand the interplay of 
structural, electronic and magnetic variables in these 
materials within an unified framework, it is technically 
daunting to 'anneal' all these degrees of freedom simul- 
taneously. We argue that it is reasonable to abstract 
the problem of structural order, solve it, and then set 
up the electronic-magnetic problem on the appropriate 
structural motif. Our paper is organised as follows. In 
the next section we present the overall model for the 
DP's, involving the structural, magnetic and electronic 
degrees of freedom, and indicate how an effective struc- 
tural model can be formally extracted. The section after, 
discusses a phenomenological model for structural order. 
We then discuss our method and computational vari- 
ables. This is followed by the results, first on 'ordering' 
for materials of the form A2BB'06, then for the 'doped' 
cases, A2Bi+a^B'i_a^06, and finally for the ternary sys- 
tems, A2BB'i_2yB"^06. We conclude with a summary of 
our results. 



II. STRUCTURAL ORDER AND MAGNETISM 

The magnetism in the DP's is intimately related to 
electron delocalisation, w hich in turn depends on the spa- 
tial pattern of B, B' ion d^ l ^ l ^ l Let us write down the 
Hamiltonian in terms of the electronic and magnetic de- 
grees of freedom to illustrate the crucial role of spatial 
B-B' order. Let B be the 'magnetic' ion and B' the non 
magnetic one. We define a binary variable 77^, such that 
rji = 1 when a site has a B ion, and rji = when a site 
has a B' ion. The ?7's will encode the atomic positions. 
In terms of these, the model for the DP's is: 
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The / operators refer to the magnetic B sites and the 
m to the non magnetic B'. and cb' are level energies, 
respectively, at the B and B' sites, A = cb — ^b' is the 
'charge transfer' energy, and t is the hopping amplitude 
between nearest neighbour B and B' ions. We have ig- 
nored orbital degeneracy in the present model. are 
the moments on the B site, J is the Hunds coupling on 
those sites, and Jaf is the antiferromagnetic (AF) su- 
perexchange coupling between B moments when two B 
ions neighbour each other. The Vat{T]} represent atomic 
interactions between the B, B' ions. 

In the simplest case of equal proportions of B and B' 
ions, and their perfect alternating arrangement, each B is 
coordinated by B' only and vice versa. The AF coupling 
does not come into play, and electron delocalisation on 
the B-B' network generally prefers a ferromagnetic spin 
configuration. This state is also highly conducting. How- 
ever, if the atomic order is imperfect and there are B ions 
neighbouring each other, two B moments get locked into 
an antiparallel configuration. This leads to a reduction in 
the overall ferromagnetic moment, and these 'antisite' re- 
gions also hinder electron transport. The magnetism and 
transport is obviously strongly dependent on the struc- 
tural motif. What decides the atomic B-B' arrangement? 
Let us look at the formal answer first. 

For large S spins, the DP model refers to electrons 
coupled to classical magnetic moments and moving in 
a background defined by {r]}. If the atomic ordering 
problem is to be isolated from this, one should 'trace 
out' the electronic and magnetic variables. The effective 
potential controlling atomic order is 

Veffiv} = -{l/P)log j VS,T^f^m}e-^" 

Computing Veff and 'updating' atomic positions accord- 
ingly is a computationally demanding task, the Monte 
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Carlo (MC) equivalent of a Car-Parinello simulation. 
There is limited information about Vat, and the 'trace' 
is technically difficult, so we construct a simple Veff{r]} 
that is consistent with the phenomenology, rather than 
attempt an elaborate ffirst principles' calculation. We 
motivate this in the next section. 



III. MODEL FOR ATOMIC ORDERING 

The ideal DP with the general formula A2BB'06 has 
ordering of B and B' atoms at the center of alternate Oq 
octahedra. Focusing only on the B and B' ions, the ideal 
structure is simply an alternation of B and B' ions along 
the X, y and z axis. However, in imperfectly annealed 
systems, there can occur antisite defect regions where this 
ordering is reversed, and two B atoms or two B' atoms 
occur adjacent to each other. The ordering of the B-B' is 
in general neither 'perfect' nor random, it is the result of 
an annealing process. While samples with high degree of 
order have been grown^ , indicating that the atomic order 
in the structural ground state should be perfect, recent 
experiments reveal an interesting trend in the degree of 
order as a function of annealing temperature. 

In the experiments by Sarma et. alJ on Sr2FeMo06 
it was observed that there was a non-monotonic depen- 
dence of the degree of Fe-Mo ordering upon the anneal- 
ing temperature, Tann, as shown in Fig.l. The different 
samples were taken from the same 'parent' material (syn- 
thesised at high temperature and quenched to a low tem- 
perature), heated to a temperature T^nn, and annealed 
there for a duration r^nn^ say. 

If there is indeed a B-B' ordering tendency in the DP's, 
the extent of order at equilibrium would be highest at 
low Tann, progressively falling off at higher Tann where 
the order is expected to be small. The downturn at low 
annealing temperature, it seems reasonable, is due to in- 
sufficient equilibriation. Our model below and results 
are based on the assumption that: (i) there is an intrin- 
sic B-B' ordering tendency in the DP's, (ii) given long 
equilibriation time, the DP's would indeed show a high 
degree of order at low T, but (iii) under typical synthesis 
conditions/annealing protocol the system only manages 
to generate correlated configurations with short range or- 
der. The annealing temperature and annealing time are 
therefore key to quantifying the structural order. 

Since the structural (dis)order seems to be 'frozen' 
at temperature, T ~ lOOOK, much above the tempera- 
ture for magnetic order (~ 400K), the qualitative issues 
in atomic ordering can be understood by ignoring the 
electronic-magnetic variables in an effective model, be- 
low. 

In the absence of detailed microscopic knowledge about 
Veff{r]}, we used a binary lattice gas model that has the 
same ordering tendency as the real materials, viz^ B-B' 
alternation, or equivalently a 'checkerboard' pattern. In 
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FIG. 1: Experimental data for ordering vs annealing temper- 
ature, obtained by Sarma et. al 



terms of the variables 77^, the simplest such model is: 

Veff{ri} = vY,m{l-rij) (2) 

with y > being a measure of the ordering tendency. 
The ground state in this model would correspond to 7^ = 
1, 0, 1, 0, .. along each axis, ie, B, B', B, B'.. Notice that 
this approach tries to incorporate the effect of complex 
interactions between the A, B, B' and O ions (as also the 
electrons) into a single parameter (See Appendix). 

This binary model is equivalent to a nearest neighbour 
Ising model, so the equilibrium physics is very well under- 
stood^. We, however, want to explore the consequences 
of different annealing protocols on this model, to examine 
the consequences of imperfect annealing. The qualitative 
ordering effects in the lattice gas model are similar in 2D 
and 3D, so we work with a 2D structure since it allows 
ease of visualisation. 



IV. METHOD 

The atomic order that emerges from Vef / can be stud- 
ied with mean field theory if the system is in equilibrium. 
However, we want to explore non equilibrium effects due 
to poor annealing, so we use a Monte Carlo technique 
to anneal the rj variables. Most of our studies involve a 
protocol where we start with a completely random B-B' 
configuration (as if quenched from very high T), and then 
anneal it at a temperature Tann, foi" a MC run of dura- 
tion Tann Moutc Carlo stcps (MCS). Since the number of 
B and B' atoms is fixed we update our configurations by 
(i) moving to some site R^, (ii) attempting an exchange 
of the atom at R^, with another randomly picked within 
a box of size centred on R^, and (iii) accepting or re- 
jecting the move based on the Metropolis algorithm. We 
have used system size L = 20, 40, and 80, update cluster 
with Lc = 4, and 10, Tann/V ranging from 0.01 - 0.3, 
and Tann ^ 200-30000 MCS. We have studied the struc- 
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ture factor, and also detailed spatial configurations, to 
quantify the extent of order. 



V. ORDERING WITH EQUAL PROPORTIONS 
OF B, B' 

This section discusses the order in systems where the 
number of B and B' ions is equal, and perfect order is 
in principle possible. In an attempt to mimic the ex- 
perimental annealing protocol, we started with random 
initial configurations at some temperature T^nn and an- 
nealed for some time Tann- The structure factor S{q) at 
the ordering wavevector, {7r,7r}, is averaged over 40 such 
initial configurations. The extent of order, quantified by 
the peak in 6'(q), is shown in Fig. 2(a) for a lattice size 
20 X 20, and in Fig.2(b), for a lattice size 40 x 40. 

There are two noteworthy features in the data in Fig. 2. 
(i) The non monotonic behaviour with Tann that we an- 
ticipated indeed shows up in the structure factor peak, 
and (ii) there is a strong size dependence of the peak in 
5'(7r,7r), varying almost by a factor of 4 between L = 20 
and L 40! 

The downturn in 6'(7r, tt) at low Tann is due to the 
inability of the system to achieve equilibrium, at short 
annealing time, when one starts with a random initial 
configuration. With increasing Tann there is of course an 
increase in the extent of order (for given Tann and L) and 
5'(7r, tt) reaches ~ 0.6 of the equilibrium value for L = 20, 
for MCS=1000. At L = 40, however, this is suppressed 
to ~ 0.1 of the equilibrium value, for the same MCS. This 
origin of this strong size dependence becomes apparent 
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FIG. 2: 71 — 71 structure factor VS Tann Oil du 20 X 20 lattiCG, 
for different Monte Carlo runlength (annealing time). 20 x 20 
lattice (top) and 40 x 40 lattice (bottom). 




FIG. 3: Non-equilibrium, Tan 
100. 
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when we examine a typical configuration at low Tann^ 
generated by a short annealing run, Tann = 100 MCS, in 
Fig.3 

The ordering is obviously imperfect, as 5'(7r,7r) sug- 
gested, more interestingly, the system actually consists 
of a few large ordered clusters with 'phase slip' between 
them. While locally these domains are well ordered, the 
bulk 5'(7r, tt) arises from the interfering contribution of 
large domains, and this 'cancellation' depends strongly 
on the system size L. In the smaller systems, L ~ 20, 
5'(7r,7r) is decided by the larger domain, and as domains 
proliferate with increasing L, there is an increasingly bet- 
ter cancellation between the out of phase domains, and 
5(7r,7r) falls. 

Since the antisite regions are at the interface of two 
ordered (but phase slipped) clusters, in what follows we 
show the domains and domain walls, rather than the de- 
tailed atomic configuration. 

In Fig. 4 one corner site in each panel is set as reference 
and the others 'coloured' in terms of their phase relation 
to it. For concreteness assume the atom at bottom left 
corner is B. Let this site be Rq = {0,0}, and index all 
sites in terms of integers, {Rix^Riy}. Then the 'correct' 
order would require that the atom at be B if Ri^ + Riy 
is even, and B' if Rix-\-Riy is odd. Formally, we just plot 
fiVi) — 77^e*^^^^^+^^^\ with a light colour for a site that 
is correctly ordered with respect to the origin, and dark 
if it is not. 

The result of this, on a 40 x 40 lattice, is shown 
in Fig. 4. The columns correspond to different Tann = 
1000, 500, 200 MCS, left to right, and the rows to 
Tann = 0.1, 1.0, 2.0, 3.0, top to bottom. At larger an- 
nealing and intermediate Tann we see a large single do- 
main dominating the configuration, while at low or high 
Tann ^ud short Tann the structurc is more fragmented. In 
the electron problem, the domains themselves are likely 
to be 'homogeneous' ferromagnetic regions, while the an- 
tisite regions, with neighbouring B-B atoms might be an- 
tiferromagnetic. 

In an attempt to quantify the size of the clusters as a 
function of Tann and Tann we studied the full structure 
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FIG. 4: Domain pattern for different Tann and equilibriation 
time. Left column, tmc = 1000, center tmc = 500, right 
TMC = 200. Rows, top to bottom, T = 0.1, 1.0, 2.0, 3.0. 
The antisite defects are at the domain boundary. 
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FIG. 6: TT — TT structure factor vs T for heating run with 
memory for different Monte Carlo steps (annealing time) 



We also tried a protocol where a system is heated 
to a temperature Tf but gradually, through a sequence 
Ti, T2, .., retaining the memory of configurations at the 
earlier temperature. The result of this slow annealing 
is similar to the earlier protocol, ie, non-monotonic for 
short Tann 5 and tends towards the monotonic equilibrium 
response with increasing Tann- For Tann > 10^, we obtain 
essentially the equilibrium profile except at very low T. 
The results are in Fig. 6. 



VI. ORDERING WITH UNEQUAL 
PROPORTIONS OF B, B' 



factor 6'(q), averaged over ~ 50 configurations, and fitted 
the ordering peak to a lorentzian of the form 

5(q ~ {tt, tt}) « A(r/7r)/((9x - irf + {qy - + T^) 

A is an overall amplitude factor and T{Tann^^ann) is a 
measure of the (inverse) width of the cluster size. This 
is shown in Fig. 5. 
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FIG. 5: The Lorentzian width in the structure factor peak, 
to estimate the domain size/correlation length. 



When the B-B' proportion is not 1:1 it is obviously 
not possible for every B to be coordinated by B' and 
vice-versa. 'Antisite' regions, rich in B or B', would ex- 
ist even at equilibrium. Our attempt, in this section is 
(i) to capture the loss of B-B' order with increasing con- 
centration of B', say, and also (ii) to study the effect of 
restricted annealing on the system. 

Our reference for (i) above are the results of D. Top- 
wal, et. al.^^ who had prepared the series of samples 
Sr2Fei+a;Moi_a^06, with proportion of Fe in excess of 
Mo. Such samples would have antisite defects even when 
the sample is well equilibriated. The XRD pattern for 
these set of samples is given in Fig. 7. They observed 
that the peak at 19.6°, related to checkerboard ordering, 
decreases uniformly for both positive as well as negative 
X. In particular, the peak is observed to be almost absent 
in the x = 0.5 sample. 

Before embarking on a MC study of this problem it is 
useful to establish the results of a simple mean field study 
of B-B' order (at equilibrium) for varying x. Fig. 8 shows 
the result of such a calculation, done via the mapping of 
the lattice gas to an Ising model at constant 'magneti- 
sation' (which corresponds to the difference in concen- 
tration of B and B'). One can capture this effect within 
mean field theory by considering an AF Ising model at 
different magnetic fields, and computing the AF order pa- 
rameter, ie, the ordering peak, and the magnetisation. 
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FIG. 7: XRD data for different Fe-Mo proportions, obtained 
by Topwal et. aP 



As seen in Fig. 8, tfie {tt, tt} ordering at low tempera- 
tures decreases uniformly with increasing magnetization. 
However, it truly vanishes only when the magnetization 
is unity, correponding to a purely B (eg. Fe) or B' (eg. 
Mo) compound. This does not correspond to the exper- 
imental situation. 




FIG. 9: Structure factor vs T (cooling runs) for different 
B-B' proportions. 



fixed, while the proportions are varied. While 5'(7r, tt) 
obviously has a nonmonotonic dependence on T^nn, the 
maximum values progressively increase with decreasing 
X as expected. In the case of annealing with memory, as 
shown in Fig. 12, the behaviour is once again nonmono- 
tonic, except that the maximum of the structure factor 
nearly always hits the equilibrium value. 
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FIG. 8: Structure factor at (tt, tt) from mean field theory for 
different B-B' proportions 

Next we did Monte Carlo simulations for the Ising 
model for different non-zero values of magnetization. The 
results for cooling runs are shown in Fig. 9. The {tt, tt} or- 
dering is found to almost vanish at x = 0.5, much before 
the mean field prediction. Interestingly, the experimen- 
tal XRD data at x = 0.5 does not seem to show the 19.6° 
peak related to ordering either. 

For any specific proportion, x = 0.1, say, the varia- 
tion of 5'(7r,7r) with T^nn is shown in Fig. 10, for various 
Tann- The Overall behaviour is similar to that for equal 
B-B' proportions, namely that the non-monotonicity de- 
creases upon increasing annealing time, and gradually 
approaches the equilibrium cooling curve. Since the equi- 
librium saturation value of the (tt, tt) structure factor in 
this case is less than 1, hence the maximum values of the 
non-equlibrium curves for any given value of MC steps is 
less than the corresponding value for the equal propor- 
tion case. 

In Fig. 11, on the other hand, the annealing time is kept 



FIG. 10: Structure factor vs T (heating run without memory) 
for X = 0.1, different Monte Carlo steps/spin. Average over 
300 initial configurations. 
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FIG. 11: Structure factor vs T (heating run without memory) 
for different x, fixed Monte Carlo steps/spin=1000. 

Fig. 13, left, shows a low temperature configuration in 
a well equilibriated sample for 2:1 proportion of B:B' 
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FIG. 12: Structure factor vs T (heating run with memory) 
for different x, fixed Monte Carlo steps/spin=1000. 



species. The ordering is disturbed by the occurence of 
large antisite patches, forced by the excess B. However, 
the species which occurs in the lesser proportion is found 
to order as much as possible, and does not form nearest 
neighbours. In this 'phase-colouring' scheme, the antisite 
regions show up as checkerboard patterns. The actual Fe- 
Mo checkerboad regions, on the other hand, shows up as 
domains of a particular colour. The black and white do- 
mains are out of phase with respect to each other, while 
they are separated by thick 'antisite walls'. 

Fig. 13, right, shows the configuration for B:B'=3:1, the 
concentration at which the ordering all but disappears. 
It is observed that the antisite patches have started per- 
colating at this concentration. 




FIG. 13: Left: Equilibrium low T (=0.5) configuration for 
large MCsteps=5000, Fe:Mo=2:l (x=0.33). Right: Equilib- 
rium low T (=0.5) configuration for large MC steps=5000, 
Fe:Mo=3:l (x=0.5). 



VII. COMPUTATIONAL CHECKS 

While the previous sections essentially summarize the 
results for the Fe-Mo ordering problem, it is necessary 
to probe the sensitivity of the results to the system size 
and the updating cluster size chosen. In this section, 
we study these systematics, thereby providing computa- 
tional indicators to the robustness of our results. In Fig 
14 and 15, we plot the dependence of domain sizes on 
system size and B-B' ratio. The three columns stand 
for L=20,40 and 80 from left to right, while the three 




FIG. 14: Domain structure for various system sizes: L — 
20, 40, 80 (left to right) and B:B ratio 1 : 1, 1 : 2, 1 : 3 
(top to bottom). This is a poorly equilibriated system, with 
Tann = 200MCS. Notice the increase in antisite disorder with 
increasing B:B' disproportionation, and the rough constancy 
of the typical domain size with increasing L. 



rows correspond to B-B' ratio 1:1,1:2 and 1:3 from top 
to bottom. It is found for all the proportions that the 
rough domain size remains the same irrespective of the 
system size, which means that there are more domains for 
larger system sizes, with more structure to them. This 
is precisely the reason for the sharp suppression of the 
nonequilibrium structure factor with increasing system 
size obseved in Section V. 

The difference between Fig. 14 and Fig. 15 lies in the 
fact that the former shows configurations obtained after 
200 Monte Carlo steps, while the latter shows configu- 
rations obtained by continuing the same run upto 5000 
MCS. It is observed that for the case of equal propor- 
tions, the domains tend to clump together and become 
more uniform as the number of MCsteps increase. For 
the case of different proportions especially Fe:Mo=2:l, i 
the fraction of antisite defects for lesser MCS is substan- 
tially larger than that for greater MCS. This is due to 
the fact that for MCS=200, there is an extra source of 
antisite defect formation in addition to the disproportion- 
ation, namely insufficient equilibriation. In other words, 
even the species which occurs in lesser proportion forms 
nearest neighbours, in contrast to the equilibrium cases 
of Fig. 13. On the other hand, for Fe:Mo=3:l, the inher- 
ent disorder due to disproportionation is so large as to 
mask out the additional nonequilibrium effects. 

In Fig. 16, the domain pattern is shown for two different 
update cluster sizes, Lc = 4 (first column) and Lc = 
10 (second column). The first row corresponds to 200 
MCS, while the second corresponds to 5000 MCS. No 
appreciable qualitative change is observed, indicating the 
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FIG. 15: The same as Fig. 14, but with much longer equihb- 
riation, Tann = 5000MCS. 



robustness with respect to this parameter, at least within 
the regimes considered. 




FIG. 16: Domain pattern with two different update cluster 
sizes: 1st column: Lc — 4, 2nd column: Lc — 10. 1st row: 
200 MCS, 2nd row: 5000 MCS. L = 40. The same initial 
configuration were used as in Fig. 14 and Fig. 15. 



VIII. ORDERING IN TERNARY B-B'-B" 
SYSTEMS 

If Sr2FeMo06 is doped with tungsten to create 
the series of compounds Sr2FeMoajWi_aj06, then one 
has a compositional/structural problem involving three 
species: Fe, W and Mo. In experiments done by 
Kobayashi et.al^ it was found that the lattice param- 
eters a and c increase with increasing W concentration, 
indicating a size mismatch. However, the more interest- 
ing observation appears to be on the ordering of Fe with 
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FIG. 17: Experimental data on B-B^ site ordering in 
Sr2FeMoa;Wi-a;06; figure taken from ReP 




FIG. 18: Structure factor vs T (cooling runs) for different 
Mo-W proportions for the Sspecies case. 



respect to Mo/W. Both Kobayashi et. al. and Sarma 
et. a/.f^ observed a larger degree of B-B' ordering in the 
heavily tungsten doped compounds. The ordering is in- 
deed observed by Kobayashi et. al. to be almost 100% 

However, no such order- 



for X < 0.6, as seen in Fig 17 



ing is observed on the B' site for the W and Mo species. 
Instead, a random mixing of W and Mo is maintained 
throughout the series. 

The above results suggest that there is an effective 
short range attractive interaction between Fe and either 
W or Mo, while there is only a much weaker interaction 
between W and Mo. However, the Fe-W attractive ten- 
dency is apparently larger than the Fe-Mo one, since the 
degree of ordering increases with increasing W concen- 
tration. Hence, we model this three-species annealing 
problem in the following manner. 

We consider a site- variable S which can take three pos- 
sible values: (Fe), 1 (W) and -1 (Mo). Then, we con- 
sider the following model hamiltonian: 

This model, a variant of the Blume-Emery- Griffiths 
model^, is a generalization of the Ising model to spin 
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FIG. 19: Structure factor vs T (heating without memory) for 
different Mo-W proportions for a fixed annealing time (2000 
MCsteps/spin). 
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FIG. 20: Structure factor vs T (heating without memory) 
for Mo:W=l:l for different anneaUng times. 



5* = 1 rather than 1/2. One can reduce this model to 
an effective Ising model in the binary limits of (Fe,Mo), 
(Fe,W) and (Mo,W), with effective exchange couplings 
given hy J -\- K — 2L^ J -\- K -\- 2L and J respectively. 
Motivated by the experimental observations, we choose 
the values of the parameters J, K and L such that the ef- 
fective energy scale for Fe-W ordering is somewhat larger 
than the Fe-Mo scale, both being substantially large com- 
pared to the Mo-W ordering scale. With such a param- 
eter choice, a Monte Carlo is performed as before, and 
the extent of B-B' site ordering (not distinguishing be- 
tween Mo and W on the B' site, is again quantified by the 
structure factor at {7r,7r}. The variation of this struc- 
ture factor with temperature for different relative Mo- 



W concentrations is plotted in Fig. [18] for cooling runs. 
The extreme W-only case of course exhibits the highest 
Tc while the Mo-only case shows the lowest one, as ex- 
pected. For intermediate Mo-W concentrations, the Tc 
for the 3-species assembly interpolates between the two 
limits. Hence, if one assumes that one has annealed the 
sample long enough to reach equilibriation at any spe- 
cific annealing temperature, the B-B' ordering increases 
uniformly with increasing W concentration, as observed 
in the experiment. However, it is interesting to note that 
the behaviour gets reversed in the off-equilibrium non- 



FIG. 21: Equilibrium low T configuration for Mo:W=l:l. 
The black squares represent the Mo, while the white and the 
yellow represent the W and Fe respectively. Notice perfect 
B-B' ordering, but no ordering on the B' site between Mo 
and W. 



monotonic case, as observed in our result for heating run 

This is because, at any 



without memory given in Fig. 19 
given low temperature, the activation probability for the 
lowest Tc compound is the highest, while that with the 
highest Tc is the lowest. Hence, the compound with the 
lowest Tc reaches equilibrium fastest. This is the reason 
for the reversal. For completenes, the heating run with- 



out memory data is also shown in Fig. 20 for different 
annealing times, keeping the Mo-W proportion fixed at 
1:1. The same transition from non-monotonic to mon- 
tonic behaviour is observed in this case also. Finally, 
Fig. 21 shows a ground state equilibrium configuration 
for Mo:W=l:l obtained from a cooling run. It is ob- 
served that there is perfect ordering on the B site, cor- 
responding to the Fe sublattice. On the B' site. Mo and 
W compete for space. There appears to be no significant 
ordering amongst these latter. Hence, overall this sim- 
ulates the experimental situation for Sr2FeMoi_a;Wa^06 
quite effectively. 



IX. CONCLUSION 

We have shown that antisite disorder can occur due to 
three reasons: (i) insufficient annealing at low tempera- 
tures leading to non-equilibrium configurations, (ii) suffi- 
cient annealing leading to equilibrium configurations but 
at high temperatures close to the order-disorder transi- 
tion, and (iii) sufficient annealing leading to equilibrium 
configurations at low temperature, but with different con- 
centrations for the B and the B' species. 
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X. APPENDIX 

We provide the mapping between the lattice gas model 
and the Blume-Emery- Griffiths model here. If we repre- 
sent a three-component lattice gas of species A,B and C 
using spin S=l, with components 1,0,-1, then the number 
of atoms of each type is given by: 

1 ^ 

i=l 
1 ^ 

^ <=1 

Nc = f2(^-Sf) (4) 

1=1 

An effective spin model can be written as: 

H = ^[JSiSj + KSfS] + L{SfSj + S]Si)] (5) 



where the three parameters J, L above are related to 
the 6 interspecies nearest neighbour interactions by: 

J = Eaa/^^Ebb/^-Eab/2 

K = Eaa/ 4 + Ebb 1 4 + Ecc + Eab/"^ — Eac — Ebc 
L = Eaa/^-Ebb/^-Eac/'^^EbcI'^ (6) 

It is to be noticed that the relation between the 'Ising' 
parameter J with the energies involve only the A and 
B species, and is identical to the corresponding relation 
for the binary alloy case. The explicit relation between 
the Ising parameter J and the lattice gas parameter V is 
given by y = 4J provided the bonds are counted only 
once (i.e., the summation over i, j is restricted to i < j). 
Hence, in terms of a lattice gas model, the three param- 
eters Eaa^ Ebb and Eab can be expressed in terms of a 
single parameter: 

V = [Eaa + Ebb — '^Eab) 
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